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ABSTRACT 

Within the next decade, ground based gravitational wave detectors are in principle capable of deter- 
mining the compact object merger rate per unit volume of the local universe to better than 20% with 
more than 30 detections. These measurements will constrain our models of stellar, binary, and star 
cluster evolution in the nearby present-day and ancient universe. We argue that the stellar models are 
sensitive to heterogeneities (in age and metallicity at least) in such a way that the predicted merger 
rates are subject to an additional 30-50% systematic errors unless these hetcrogenities are taken into 
account. Without adding new electromagnetic constraints on massive binary evolution or relying on 
more information from each merger (e.g., binary masses and spins), as few as the ~ 5 merger detec- 
tions could exhaust the information available in a naive comparison to merger rate predictions. As 
a concrete example immediately relevant to analysis of initial and enhanced LIGO results, we use a 
nearby-universe catalog to demonstrate that no one tracer of stellar content can be consistently used 
to constrain merger rates without introducing a systematic error of order 0(30%) at 90% confidence 
(depending on the type of binary involved). For example, though binary black holes typically take 
many Gyr to merge, binary neutron stars often merge rapidly; different tracers of stellar content are 
required for these two types. More generally, we argue that theoretical binary evolution can depend 
sufficiently sensitively on star-forming conditions - even assuming no uncertainty in binary evolution 
model - that the distribution of star forming conditions must be incorporated to reduce the system- 
atic error in merger rate predictions below roughly 40%. We emphasize that the degree of sensitivity 
to star-forming conditions depends on the binary evolution model and on the amount of relevant 
variation in star-forming conditions. For example, if after further comparison with electromagnetic 
and gravitational wave observations future population synthesis models suggest all BH-BH binary 
mergers occur promptly and therefore are associated with well-studied present-day star formation, 
the associated composition-related systematic uncertainty could be lower than the pessimistic value 
quoted above. Further, as gravitational wave detectors will make available many properties of each 
merger - binary component masses, spins, and even short GRB associations and host galaxies could 
be available - many detections can still be exploited to create high-precision constraints on binary 
compact object formation models. 
Subject headings: Stars: Binaries: Close 

1. INTRODUCTION 

Ground based gravitational wave de- 
tector networks (LIGO, described in 

lAbbott et al. (The LIGO Scientific CoUaborationll 
(|2003f l: and VIRGO, at the Virgo project web- 
site www.virgo.infn.it) are analyzing the re- 
sults of a design-sensitivity search for the sig- 
nals expected from the inspiral and merger of 
double compact binaries (here, NS-NS, BH-NS, 
BH-BH) (extending, for example, the search in 
Abbott et al. (The LIGO Scientific Collaboration)! 
20081 ). Sensitivity improvements in LIGO and other 
interferometers are expected over the next decade that 
will make multiple detections a near certainty. For 
example, based on the short lifetime of the very massive 
black hole X-ray binary IC-10 X-1, iBulik et all (pOOSh 
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predict that even the current generation of interferome- 
ter has a good chance of detecting a (high-mass) merger. 
With moderate improvements in detector sensitivity 
that will be in place by mid-2009 ("enhanced LIGO"), 
multiple detections are plausible. More conservatively, 
theoretical calculations which e xplore a wide range 
of stil l-plausible assu mptions [O'Shaughnessy et all 
(|2008bD (PS-GRB) and lO'Shauehnessv et all ^2008^ 
(PS-E)] predict that the advanced LIGO network is 
likely to detect several tens of mergers per year, allowing 
the merger rate per unit volume to be determined in 
principle to within 20%. In fact, advanced LIGO can 
determine the merger rate p er unit volume significantly 
more precisely (20%; see.e.g. lO'Shaughnessv et aI]l2007L 
(PS-E2)) than measurements have constrained the star 
formation history of the Milky Wa y and the distant uni- 
verse ( often at least 3 0%; s ee.e.g. iPerez-Gonzalez et al.l 
120061 IWilkins et all l2008l and references therein). 
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against which the predictions of lO'Shaughnessv et al.l 
((20083) and other theoretical models are normalized. 
As theoretical predictions can be no more precise than 
their input, even though a large number of merger 
detections are likely, advanced LIGO measurements 
cannot distinguish between different hypotheses about 
how merging binaries are produced if those merger rates 
differ by less than 0(30%), on the basis of the number 
of mergers alone. 

With the ability to measure both the number and 
properties of merging compact binaries, LIGO has long 
been expected to provide invaluable assistance in better- 
constraining hyp otheses regarding corn pact binary for- 
mation (see, e.g., ISadowski et al.l [20081 and references 
therein). Given a systematic error of e in any merger 
rate prediction, only approximately unique detec- 
tions are needed to determine if reality is consistent with 
a model. Even optimistically assuming that three types 
of binaries can be distinguished (e.g., from their com- 
ponent masses) and their rates estimated at this level 
of accuracy, the fraction of a priori plausible models still 
consistent with the three merger rates provided by LIGO 
is comparatively large: ~ [(logl -I- e)/2]'^ (assuming two 
orders of magnitude uncertainty in the a priori plausi- 
ble merger rate for each type of binary), at best com- 
parable with the minisculc fraction of parameter space 
needed to constrain a high-dimensional model weakly. 
For example, for the seven-dimensional binary evolution 
models compared with obser vations of double pulsars in 
lO'Shaughnessv et al.l ()2008c| ). each parameter would be 
constrained to just [(logl -I- e)/2]^/'' ~ 0.3 of its a pri- 
ori range using merger rate estimates alone, comparable 
(albeit complementary) to the information provided by 
electromagnetic observations of double pulsars. On the 
contrary, had systematic errors been smaller, then detec- 
tion of n binaries of each of 3 types should imply an ac- 
curacy (21nl0)-3/7n-3/i4 0.2(n/100)~3/i4 in each pa- 
rameter. Furthermore, since this systematic uncertainty 
is introduced through our lack of knowledge about the 
nearby and ancient universe, even though third genera- 
tion detectors such as the Einstein Telescope will harvest 
vastly more mergers, they will be similarly limited when 
comparing their observed merger rates with theoretical 
models that rely upon existing surveys of star forma- 
tion. To take full advantage of the many mergers that 
in-construction and third-generation instruments will de- 
tect, compact-object theorists will need to compare the 
distributions of binary parameters expected from theory 
(i.e., masses, spins) with observations. 

In this paper we estimate the limiting systematic er- 
ror introduced into any theoretical prediction of binary 
compact object merger rates through the star formation 
history of the universe. We furthermore explain that the 
relevant uncertainty is not merely overall normalization 
of the nearby and even distant star formation history. 
Instead, we argue that merger rates, particularly binary 
black hole merger rates, can also be sensitive to the corre- 
lated distribution of age and metallicity of their progeni- 
tor star-forming regions. Though high-precision surveys 
and spectral energy distribution (SED) reconstructions 
of galaxies may precisely determine the mean star forma- 
tion rate and metallicity by the epoch of advanced LIGO, 
(;see.e.g. lRenzinill2006HHoDkins et al.ll2003l : IWilkins et all 



l2008l and references therein) , the more delicate analyses 
which estimate the distribution of star forming condi- 
tions, particularly those of low metallicity that are far 
more apt to produce ma ssive black hole bin aries, remain 
in their infancy (see.e.g. jPanter et al.ll2008l ). Given that 
advanced LIGO and future gravitational-wave detectors 
could observe mergers produced from binary stellar evo- 
lution out to as far as z ~ 2 (e.g., for an optimally ori- 
ented 30 -I- 3OM0 binary black hole merger), an epoch 
of rapid star formation in massive galaxies, the relevant 
composition distribution needed to eliminate this sys- 
tematic uncertainty is unlikely to be available in the near 
future. Equivalently, gravitational-wave interferometers 
will soon provide a uniquely accurate and potentially 
uniquely biased probe into the formation and evolution of 
high-mass stars in the early and low-metallicity universe. 

1.1. Outline and relation to prior work 

As discussed in more detail in ij2l to account for 
local-universe inhomogeneities and to simplify the 
intrinsically mass-dependent results of gravitational 
wave searches, previous searches for gravitational wave 
inspiral and merger waveforms have "normalized" their 
result by the amount of blue light within the relevant 
t ime- averaged detection volume; see the discussion 
in lAbbott et al. (The LIGO Scientific Collaboration)! 
(j2008f ) as w ell as the consider a bly in ore detailed pre- 
sentation in iFairhurst fc Bradvl l)2007D . By choosing to 
express results as a "merger rate per unit blue light," 
however, the authors limit the accuracy of any attempt 
to compare merger rate predictions with their observa- 
tions: as emphasized above, such a comparison is helpful 
only to the level of accuracy that "mergers per unit blue 
light" can be uniquely defined. The systematic error so 
introduced is unlikely to seriously limit astrophysical 
comparisons once detections are available in the near 
future: initial and enhanced LIGO results, expected to 
have at best a handful of detections, will not reach this 
level of accuracy. But this composition-based systematic 
error is comparable to several formal uncertainties often 
quoted in relation to already-published upper limits and 
is therefore already relevant to anyone attempting to 
constrain their models with existing observational upper 
limits. In short, anyone planning on using or expressing 
results in this form should be aware of its limitations. 

That being said, at present, gravitational-wave detec- 
tors survey only the nearby universe, where uncertainties 
in the distances to galaxies dominate over photometric 
error s (74% vs 31%, respectively; see iKopparapu et al.l 
|2008[ ). At best, the cumulative asymptotic luminosity 
can be determined only to 0(10%). All of these uncer- 
tainties arc comparable or greater than the uncertainty 
introduced into any astrophysical interpretation by as- 
suming the number of mergers is proportional to blue 
light. The uncertainties discussed in this paper there- 
fore bound below the accuracy of any comparison between 
merger rate predictions and observations. 

These limiting uncertainties arise because most sim- 
ple prediction (or "normalization") methods build in an 
implicit assumption of homogeneity of star-forming con- 
ditions. But because binary mergers are rare and ex- 
ceptional events themselves and naturally arise more fre- 
quently from rare and exceptional conditions (e.g., old 
star formation or low metallicity), assuming homogene- 
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ity builds in systematic errors greater than the hmiting 
uncertainty desired for advanced detectors. As outhned 
above and described in § [21 we argue that dividing the 
rate of mergers by the amount of blue light in the de- 
tection volume oversimplifies the (implicit) inverse prob- 
lem: predicting how many mergers should occur given an 
amount of blue light. To demonstrate that other bands 
give different yet potentially equally relevant normaliza- 
tions, we introduce a multi-band galaxy catalog for the 
local universe. In ^ we demonstrate that, after a star- 
burst, different models of binary evolution and different 
types of binaries lead to different conclusions about the 
time-dependent ratio of mergers and light. We use this 
tool to estimate the systematic error introduced by nor- 
malizing to blue light or, more generally, any single-band 
normalization. For advanced detectors, galaxy catalogs 
will not be available. Nonetheless, as demonstrated in §2] 
the gravitational wave detection rate should not be cav- 
alierly normalized to the mean properties of the universe 
on large scales: exceptional circumstances (here, metal- 
licity) can introduce systematic errors at least as large 
as the limiting uncertainty expected of advanced detec- 
tors. Though the exact magnitude of the effect cannot 
be determined without an equally exact theory of bi- 
nary evolution, wc estimate that even modestly reliable 
predictions could require fairly detailed input regarding 
the composition of the universe within reach. Finally, 
to clearly illustrate the effects summarized in this paper, 
in Sj5] we show that concrete, plausible predictions for 
a two-component universe cannot be well-modeled by a 
time-independent or homogeneous one. 

2. LOCAL GALAXIES IN MULTIPLE BANDS 

In the past, the number of mergers per unit blue light 
has been used to normali ze the sensitivity of searches 
and interpret upper limits (jPhinnev 19911. Because blue 
light roughly traces current star formation and because 
many double neutron star mergers occur fairly soon after 
their progenitor binary's birth, this ratio was expected to 
be proportional to the fraction of massive stars that, af- 
ter their second supernova, form bound double neutron 
stars that merge within a Hubble time. This assump- 
tion was applied wide ly in the theoretical (lPhinnevlll991l: 
iKalogera et al.|[200ll) and expe rimental (jNutzman et al.l 
(|2004D : [Kopparapu et all ()2008D henceforth LGC) litera- 
ture. The blue light density locally and at moderate red- 
shift can be measured very accurately (30%, dominated 
by per-galaxy distance errors). Being larger than the 
detector's intrinsic systematic error target (15%), this 
measurement error has implicitly been treated in the 
gravitational-wave literature as the relevant systematic 
uncertainty on binary merger upper limits per unit star 
forming matter (LGC). 

While adequate to zeroth logarithmic order, the tra- 
ditional approach is accurate only to the degree that 
the universe satisfies two approximations: (i) that only 
present-day star formation dominates the present-day 
compact-object coalescence rate and (ii) that all galaxies 
are sufficiently similar that twice as much blue light cor- 
relates directly with twice as many mergers. In reality 
elliptical galaxies are expected to contribute a significant 
proportion of all present-day compact bina ry coalescence 
detections ()de Freitas Pacheco et al.ll2006[ ). particularly 
from BH-BH binaries (PS-E). Because elliptical galax- 



ies formed their stars long ago and under different star- 
forming conditions than the stars which produce most of 
the present-day blue light, normalizing coalescence de- 
tections to blue light misrepresents the relevant degrees 
of freedom and loses information. 

To provide concrete scenarios to demonstrate that old 
star formation and differences between galaxies (e.g., 
between ellipticals and spirals; between galaxies of dif- 
ferent formation history and metal content) can signif- 
icantly influence present-day star formation and to as- 
sist in re-evaluating the systematic error associated with 
present-day short-range gravitational-wave observations, 
we construct a local-universe galaxy catalog with more 
information than just blue light. Rather than use a 
deep survey with limited sky coverage to investigate 
the properties of galaxies in the large distance limit, 
given the range relevant to the current generation of 
gravitational-wave interferometers we choose to extend 
the previous B-band catalog provided in LGC. Ideally, 
we would demonstrate the importance of both inhomo- 
geneity and old star formation by using spectra of all 
relevant galaxies (e.g., all within ~ 160 Mpc for BH- 
BH mergers for initial LIGO) to reconstruct their star 
formation and composition histories and convolve each 
with an appropriate model for binary evolution. Though 
the situation may change as sky coverage of large-scale 
surveys improve, at present only photometric informa- 
tion is available for all galaxies out to the Virgo cluster. 
We have used the LEDA catalog to extract corrected 
U, V, and B apparent magnitudes, best distance esti- 
mates, and morphological classifications for ~ 38, 000 
galaxies; we convert these magnitudes to luminosities 
using the zero-point conventions adopted in the Ap- 
pendix. Though an extensive literature exists address- 
ing methods with which to reconstruct star formation 
histories, metallicities, and extinct ions from photometric 
and spectr al obs ervations (see,e.g..lGallagher et al.l | 1984 ; 
Kennicuttlll998tlHopkins et al.ll2003l:ICalzetti et al.boOT : 
Perez-Gonzalez et al.ll200d and references therein) , with 
so few bands we cannot reliably invert and reconstruct 
detailed properties of our galaxy set, even assuming the 
catalog uses a good IR correction to reconstruct the in- 
trinsic U, B,V magnitude from highly obscured star for- 
mation. At best we would be limited to an 0(30%) sys- 
tematic error in the star-fo rmation history recons truc- 
tion (see, e.g., Table 4 in Pere^ Gonzalez et aIll2006D . We 
therefore work directly with the published corrected lu- 
minosities. By way of example, Figure[l]shows the cumu- 
lative luminosity versus distance for three of the bands 
provided in the catalog. At large distances, these three 
quantities match onto the average values per unit vol- 
ume estimated from local-universe cosmological surveys, 
as discussed in the Appendix. 

Ignoring differences in when and how galaxies form 
their stars introduces a systematic error which can be 
simply (under)estimated by comparing the fraction of 
stellar mass and blue light due to all morphological el- 
liptical galaxies inside our detection volume (Figure [1]). 
At large distances, elliptical galaxies account for 60% 
of all stellar mass but 20-40% of all light, depending 
on the band used; this well-known difference is exten- 
sively describe d in the h istorical and pedagogical liter- 
ature fsee.e.g.. ICoxll2000l ). This band-dependent differ- 
ence immediately implies that any merger rate predic- 
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Fig. 1. — Top panel: Cumulative luminosity versus distance for 
three different bands (U,B,V, shown as gray, blue, and purple), 
shown per unit Lio = IO^^Lq. Also shown are (i) a thin blue 
lino corresponding to the large-distance limit predicted from back- 
ground light (sec the Appendix) and (ii) a thin vertical black line 
at the approximate completeness limit of the survey. The two more 
rapid steps in cumulative luminosity correspond to the local group 
(at a few Mpc) and Virgo cluster (at 20 Mpc). Because our cat- 
alog does not have U,V-band light for all galaxies in our B-band 
complete catalog, the large-distance behavior of these bands docs 
not approach L oc at large distances. Bottom panel: Fraction 
of all U,B,V light inside d contributed by elliptical galaxies, versus 
distance; the same color scheme is used. 



tion Rpredict bascd on multiplying the total amount of 
light times some merger rate per unit light must have a 
systematic uncertainty of order this composition uncer- 
tainty: 
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■predict — f light R-sp ^" (1 



flight)Rel 
Rtrue — fRsp + (1 ^ .f)Rel 

5R={Iught - f){Rsp - Rei) ^ 6fO{R) (1) 

where Rsp^Rei are merger rates per unit star- forming 
mass of spiral and elliptical galaxies respectively, / is 
the mass fraction in spirals, and fught is some (band- 
dependent) light fraction in spirals. Assuming composi- 
tional or age differences cause one or the other population 
to dominate the present-day merger rate, the system- 
atic uncertainty introduced by assuming the light con- 
tent traces mergers should be at least 5f = 0.6 — 0.4 = 
20%. As discussed in the appendix, similar uncertainties 
are obtained if reasonable a priori mass-to-light ratios 
are adopted for the two morphological types (see, e.g., 
iLipunov et al.l [T9950 or if more sophisticated estimates 
for M/L are adopted rather than a simple morphologi- 
cal cl assification (see, e.g.. Figure 18 in lBlanton fc Roweii 
l2007t) . To get a better estimate and to determine what 



normalization is relevant - mass, light, or some combi- 
nation thereof - we must model the relative proportion 
that past and present star formation in elliptical and spi- 
ral galaxies produce present-day mergers. 

3. LIGHT AND MERGERS LAG STAR FORMATION 

Since gravitational radiation drives merging binaries 
together exceedingly slowly, particularly for binaries with 
black holes which arc likely not kicked close together in 
supernova explosions, binaries born many Gyr ago in 
now-old stellar populations produce a significant and oc- 
casional ly overwhelming fraction of all pr esent-day merg- 
ers (see Ide Frei tas Pachcco et alj (|2006l), a s well as the 
discussion in [O'Shaughnes sv et alTl^OOSbf ) and PS-E). 
The ratio of mergers to light in that galaxy will there- 
fore depend not only on the star formation history of the 
galaxy but also on the relative rate of decay of mergers 
and light after a burst of star formation. The latter is sig- 
nificantly model- and binary type-dependent; see Figure 
[3] as well as the more detailed examples in PS-GRB and 
PS-E. No one normalization will work perfectly for all 
assumptions about binary evolution; for example, blue 
light and binary black hole mergers will rarely evolve at 
the same rate. 

3.1. Estimating systematic error for blue light 

The systematic error introduced by choosing to nor- 
malize to blue light instead of a quantity that decays as 
the desired (model-dependent) merger rate can be esti- 
mated by monte carlo over a large array of binary evo- 
lution simulations and a range of galaxy star formation 
histories. Specifically, if for simplicity we assume all star 
formation occurs in similar conditions, flux in various 
bands fu,g,... as well as the present-day total mass and 
star formation rate dp/dt can all be expressed as a con- 
volution: 

jx{0)^ldTKx{0-T)^{T) (2a) 
Rnit)^ f dTKD{t~T)^{T) (2b) 



where jx is the luminosity density emitted per unit vol- 
ume in band X and Rjj is the detection rate for a network 
of gravitational- wave detectors o f some fixed sensitivity; 
see lO'Shaughnessv et al.l ()2008af l for details. The kernels 
Kx can be extracted from siniple st ellar population li- 
braries fsee lBruzual fc Charlotll2003l for details and the 
appendix for a summary); we obtain the kernels Kjj from 
PS-E. To incorporate the influence of old star formation 
on both merger rates and present-day galactic luminosi- 
ties, we explored a one-parameter model motivated by 
studies of the star formation history of the universe: 
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where e is a dimensionless parameter indicating the rel- 
ative importance of old star formation, t = is the 
present, t ~ — T is the big bang, T = 13.5 Gyr is the 
age of the universe, Tq = 1.5 Gyr is a characteristic de- 
cay time chosen so that the shape of the star formation 
rate (SFR) reproduces the large-redshift peak in the cos- 
mological SFR at large e (and in general resembles the 
overall cosmological SFR at, e w 30; see Figure [2]), and 



Po = lA/0yr~^ is a characteristic value for a galaxy's 
star formation rate. 

Given the SFR model [Eq. [3], which depends linearly 
on e, and the definitions of Eq. [2l all these quantities 
X{e) depend linearly on e. The slopes dX / de/ X{{)) = 
mx tell us the relative importance of old versus young 
star formation to the quantity X] a larger m implies 
greater sensitivity to old star formation. To give a sense 
of scale, because the star formation history of the uni- 
verse resembles p(e ~ 30), values of m > 1/30 ~ 0.03 im- 
ply very strong dependence on old star formation: with- 
out old star formation, the quantity X would be at least 
a factor 2 smaller. To use a very concrete example, lu- 
minosities scale as 

Lx{e)=Lxm^ + mxe) (4) 
for some mx that we can calculate by evaluating Lx (i.e., 
by convolution with Kx) for any two unequal e. Simi- 
larly, for each population synthesis model and each type 
of binary q, we can calculate niq (e.g, q =BH-BH,BH- 
NS,NS-NS). The distribution of (say, msH-BH ■, 
which for clarity we will denote by ttibh) then indicates 
the range of sensitivities that q binaries can have to old 
star formation. 

Figure [3] shows our results for the distribution of m, 
both for the various types of light (vertical bars) and 
mergers (distributions, sampling a range of binary evo- 
lution assumptions). A specific set of star-forming con- 
ditions (preferred values for m), time-evolution history 
(preferred e), and mass completely characterize that 
galaxy's present-day observables. To be concrete, that 
galaxy contributes to the cumulative blue light and num- 
ber of BH-BH detections as 

alaxy 

(1 + emBH) 

Lb:Lb{0) galaxy 

(1 + emB) 

where the leading-order term is proportional to the mass. 
Because we assume all star-forming conditions arc simi- 
lar, the values of ni are the same for all galaxies and the 
cumulative detection rate and light inside a volume can 
be found by summing over all: 

RD = RD{0){l + {e) tubh) (5a) 
LB=LBmi + {e) tub) (5b) 
where (e) denotes the mass-weighted average e needed 
to reproduce the globl star formation history of the uni- 
verse and thus where for simplicity we further assume 
(incorrectly) that a galaxy's star formation history is in- 
dependent of its mass. Therefore the ratio M = Rd/Lb 
of BH detections to blue light will explicitly depend on 
the model-dependent factor ttibh- 

When blue light and mergers have exactly the same 
delayed response to star formation, the ratio Rb/Lb 
is totally independent of the star-formation history and 
therefore provides an excellent tool with which to con- 
strain the underlying theory of binary evolution. In our 
notation, when vriB = itibh the two factors in Eq. [5] 
cancel, leading to a ratio that is independent of e. More 
generally blue light and mergers do not mirror one an- 
other. Adopting a blue light normalization Nb by as- 
suming niBH-BH — * iTT-B in Eq. [S] introduces a bias. To 
be explicit, blue light normalization assumes the once- 
and-for-all proportionality 
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Fig. 2.— A plot of {dp(t,t)/dt)/{dp{T)/dt) (Eq. [Sj, the one- 
parameter star formation history model adopted in the text, for 
e = 0, 10, 30 (blue, yellow, re d, respectively). Also sh own are mod- 
els f or p(f)/p(0) drawn from INagamine et al.l I I2006I ) (red, dotted) 
and ISpringel fc Hernguistl 1)20031 ) f black, dotted). Near e = 30 
our one-parameter model reasonably mimics the time dependence 
of the st ar-formation history of the universe as well as of massive 
galaxies IIHeavens et al.l2004l) : near the present, the model is nearly 
e independent. The sensitivity of predictions such as Eq[2]to e, as 
measured by "slopes" m, tell us about the relative impact of old 
versus young star formation. 



A more detailed model that allows blue light and mergers 
to have different delay kernels K has a different normal- 
ization factor ; mBn), which is greater than Nb by 
a bias factor 



bias 



A/'((e) ; niBH) _ (1 + (e) mBn) 



Nb 



(1 + (e) tjib) 



(7) 



This bias varies depending on the model being studied. 
Figure [3] implies that the most-likely values for fbias are 
between 1.4 (BH-BH) to 1.2 (NS-NS) based on a pre- 
ferred value (e) = 30 mentioned above and in Figure ^ 

Bias isn't the most pertinent problem, however; we can 
always eliminate it by adopting a different convention for 
Rd/Lb that corresponds to the results predicted by a 
"typical" model. To continue with the example above, we 
can adopt a "typical" normalization Nav corresponding 
to Eq. ([U but with mBH {rriBH)- By using such a 
typical model, the relative bias gbias between Nav and TV 
can be much reduced: 



sifTlBH-BH) — 



(1 



= 1 



(1 + (e) {mBH) 

(e) {mBH - {mBH)) 



(1 + (e) {mBH) 



(8) 



Nonetheless, even if we adopt the single best ratio for 
Rd/Lb, fluctuations between models are still sufficiently 
significant to significantly influence results. Specifically, 
the variance a of In gtias is 



1 



(mBH) e 



(9) 



which based on Figure [3] can be of order 1.26 (BH-BH, 
NS-NS) to 1.17 (NS-NS) at one standard deviation. We 
conclude that, depending on the type of binary involved, 
comparisons between theoretical models and any single, 
model-independent quantity i?£)/LB inevitably introduce 
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Fig. 3. — Systematic error due to optical merger tracers: Top 
panel shows the distribution of m for BH-BH (red), BH-NS (green) 
and NS-NS (blue) mergers, as well as the m values predicted for 
U (gray), B (blue), and V band (black) light. Note that since our 
simulations contain many NS-NS mergers that occur soon after 
formation, U-band light provides the most reliable tracer for NS- 
NS merger rates. Bottom panel shows the distribution of relative 
systematic error x = log 5J\f/Af introduced by normalizing to blue 
light, as predicted from the distribution of m using Eq. [7] No one 
band can reproduce all merger rates for all of the one-parameter 
star formation histories; typical systematic errors introduced by an 
inappropriate normalization are 0(20%). 



a > 30 — 40% systematic error into comparisons with 
binary evolution models at 90% confidence. 
Is this bias really a problem?: The above calculation 
seems to suggest that, given a binary evolution model, 
blue light normalization of merger rates is biased by a 
known and easily-calculable factor [Eq. [?]• This cor- 
rection factor can be calculated and removed post-facto, 
when rate predictions are compared with observations. 
In other words, no bias need be introduced by normaliz- 
ing to any mass or light measure, so long as we can con- 
fidently relate that measure to the present-day merger 
rate, given assumptions about how binary mergers lag 
star formation of different types. 

Additionally, at large distances the universe becomes 
homogeneous; all different light tracers become propor- 
tional, removing the need for choosing a preferred mass 
tracer. Normalization is most naturally made per unit 
volume; rate predictions are made on the basis not of 
galaxy models but on the star formation history of the 
universe (PS-GRB, PS-E and references therein). In this 
asymptotic case normalization is apparently unambigu- 
ous and model-dependent corrections can be reinserted 
later. 

In fact, as we show below, independently of delay time 



corrections, fluctuations in composition also introduce at 
least as significant an uncertainty. The elliptical galaxies 
that host the most extreme metallicities are known to 
form their stars extremely early. While we could indeed 
correct for the contribution of old stars if all star forming 
conditions were similar and if the star formation history 
of the universe was sufficiently well-known, in the real- 
istic heterogeneous universe systematic uncertainties in 
delay time and composition must both be included. 

4. HETEROGENEITY AND BIAS 

Star-forming conditions are known to be highly hetero- 
geneous in time as stars gradually process metals within 
a galaxy, particularly for less m assive galaxies whic h 
undergo extended star formation (jHeavens et al.|[2004D . 
Even at the present e poch star forrning conditions vary 
dramatically (see.e.g.. iGallazzi et al.ll2008HPanter et al.l 
120081 and reference s therein). Both lPanter et all (|2008f) 
and TGallazzi et al] ()2008D have concluded (in their Fig- 
ures 6-8 and Table 6, respectively) that nearby galax- 
ies are likely to have metallicities Z with log Z/Zq be- 
tween —0.5 and 0.2. Young star-forming galaxies have an 
even broader range of metallicitie s, with logZ/Zp) be - 
tween -1.5 and 0.2 (Figure 10 of IGallazzi et al.ll2008f l. 
Though some authors have suggested even more signif- 
icant differences, such as a tendency towards producing 
mor e massive stars than usua l (a "top heavy IMF" ; see, 
e.g. llo^ ms fc Beacom|[200l for a discussion of models 
and observational constraints), and though an increased 
number of massive stars should correspondingly increase 
the detection rate of compact binary coalescences, in this 
paper we conservatively limit attention to the more well- 
constrained issue of metallicity fluctuations. 

The gravitational- wave detection rate Rd depends sen- 
sitively on the metallicity of the gas from which the 
progenitor binary stars form, as metallicity influences 
their structure and binary evolution. For example, 
observations of massive stars have demonstrated that, 
as expected given the larger photon cross-sections of 
metals over hydrogen, massive stellar winds in crease 
signifi c antly with more metal c ontent (see, e.g., IVinkI 
(|2008f) . ISchr6der fc Cuntl l|2005f) and references therein). 
Wind loss determines the relation between initial stel- 
lar mass and final compact remnant mass of indi vid- 
ual stars (see, e.g., Figure 1 in iBelczvnski et 311120021 ): as 
both the likelihood of a progenitor of mass and the 
volume inside which a compact binary of chirp mass 
M can be observed depend sensitively on mass, metal- 
licity fluctuations are expected to lead to significant 
changes in the relative likelihood and detectability of 
compact binary mergers. Metallicity could also influ- 
ence binary evolution in other ways, such as the amount 
of mass lost during nonconservative mass transfer or a 
common-envelope phase. Unfortunately, neither obser- 
vations nor theory provide an unambiguous answer for 
the magnitude of the effect. Theoretical methods rely 
on many unknown phenomenological parameters to char- 
acterize complex physical processes such as common- 
envelope evolution. Not only do these many unknown 
parameters influence merger and detection rates by or- 
ders of magnitude (Belczv nski et al.l l2002| ) , they do so 
in a highly-correlated fash ion (see, e.g., Appendix B in 
lO'Shaughnessv et al.ll2005[ ). Generally speaking no sin- 
gle parameter, including metallicity, produces an unam- 
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biguous trend everywhere in the parameter space. And 
equally generally the trends relevant for one type of bi- 
nary (BH-BH, say) often bear little relation to the trends 
for other types, particularly after marginalizing over one 
or more other parameters. 

Despite these challenges, we can fairly easily estimate 
the order of magnitude of the systematic error intro- 
duced by ignoring heterogeneity. As a first approxima- 
tion we assume the composition of the universe is time- 
independent and estimate the present-day merger rate, 
averaging over the heterogeneous local universe's metal- 
licity distribution p(logZ), as 



{Rd)^ d\ogZp{\ogZ) 



dN 
dtdVr 



X / dMp{M\Z)Vc{M) 



dlogZ p{\ogZ)RD{Z) 



(10) 



(11) 



where log Z is the log of the metallicity; dN/dtdVc is the 
merger rate in these conditions due to all past star forma- 
tion (and implicitly includes an integral over all time); 
p(log Z) is the fraction of star formation occurring in 
those conditions; p{M\Z) is the (chirp) mass distribution 
of merging binaries formed due to Z\ and Vc{M) is the 
detection volume for binaries of (chirp) mass M, which 
we estimate using the usual power-law formula and an 
estimate Cy of the range at which a gravitational-wave 
network can detect a s ingle d ouble neutron star inspiral. 
lO'Shaughnessv et all ()2008bf ) and lO'Shaughnessv et all 
(l2008aD have previously performed calculations of Rd{Z) 
for a range of metallicities and binary evolution assump- 
tions. Based on their raw data, we estimate that the 
primary trend due to metallicity can be characterized by 
a single first-order parameter 5 



log Rd{Z)^ log Rd{Zq)+5 log Z/Zc 



(12) 



defined individually for each type of binary and which 
allows for both the change in merger rate and in charac- 
teristic mass with metallicity. Adopting this parameter, 
the relative error made by ignoring heterogeneity should 
be of order the average value of a power law Z: 



{Rd)^{{ZIZqY)Rd{Zq) 



(13) 



iYiyZj rnax j ^rain ) 

where in the second line we assume log Z is uni- 
formly distributed between a lower and upper bound 
and conservatively adopt log Z^ax/ZQ = 0.2 and 
log Zjnin/ZQ = —0.5. Unless simulations lead to a re- 
markably metallicity-independent detection rate (i.e., S 
is very close to zero) , this expression implies that hetero- 
geneity introduces a systematic error of order 30% — 60% 
for S € [—3,6].^ Though this relative change is ex- 
tremely small compared to the differences between cur- 
rently plausible binary evolution models for merger rates, 

^ Though our calculation suggests that when <5 < the system- 
atic error would be a factor 2, when we ad opt a gaussian metalli city 
distribution which reproduces Table 6 of lGallazzi et all l|2008l ) we 
estimate a systematic error within the range stated. 



and though this uncertainty may even be smaller than 
the difference between our best StarTrack model and re- 
ality, this error is significantly greater than the target sys- 
tematic error of the LIGO analysis and greater than the 
eventual uncertainty of advanced LIGO measurements. 
What is S?: A worst-case estimate can be 
quickly extracted from the figures and results of 
lO'Shaughnessv et al.1 (|2008bf ). Merger and detection 
rates due to "elliptical" galaxies, in which the metallicity 
was varied, changed by 3 orders of magnitude (95% 
confidence). Assuming all this change was produced 
only by metallicity variation and noting metallicity 
varied by 0.5 in logZ, we expect \6\ < 6. In reality much 
of the observed variation is due to other parameters such 
as supernova kicks which strongly influence the merger 
rate. For example, a set of BH-NS merger rate estimates 
in which only Z differed suggests (5bh-ns — —2. 

Unfortunately our calculations also suggest that the 
derivative dlogRn/dlogZ changes depending on the bi- 
nary evolution assumptions adopted. For this reason, 
until a model of binary evolution can be uniquely de- 
termined, the resulting heterogeneity-dependent effect is 
at best an unknown systematic error rather than a cor- 
rectable bias. For this reason, we limit ourselves to the 
above estimate of order of magnitude of the error intro- 
duced by omitting heterogeneity in detection rate esti- 
mates. Future investigation could very well demonstrate 
that binary evolution is much less sensitive to metallic- 
ity than the above estimate; under these circumstances, 
the error introduced by ignoring heterogeneity would be 
much reduced. 

Using strong Milky-way constraints to eliminating het- 
erogeneity bias?: Observations of Milky Way compact 
binaries have long been used as stringent tests of bi- 
nary evolution. For example, attempts to explain 
the existence of individual d o uble white dwar f s (see , 
e.g.. iNelemans fc ToutI (l2005l) . iD'Antona et~al] (|2006[) . 
Ivan der Sluvs et all (l2006fl . and r eferences there i n), b i- 
narv pulsars (jWiiers et"aLl (|1992[ ). [Willems et all (|2006l) . 
and re ferences there in) , a nd X-ray binaries (see e.g. 
iPodsiad lowski et al] (|2002D a s weU as the articles and 
references in iLewin et al.l (|1995[ )) have constrained 
common-envelope evolution and the strength of super- 
nova kicks. Similarly, the challenges of reconciling 
the theoretical and observed statistics of compact bi- 
na ry populations ( c ompa re, for example, _ 
or iBelczvnski et al.l (I2002D with iKalogera et 
have also suggested constraints (jO'Shaughnessv et al.1 
120083 . henceforth denoted PSC2). 

Conceivably such strong constraints could uniquely 
determine the binary evolution model appropriate to 
the Milky Way. Combined with an understanding of 
metallicity-dependent single star evolution, we can imag- 
ine uniquely determining Rd{Z). Therefore, in an ideal 
world, by combining Rd[Z) with the metallicity distribu- 
tion of the time-evolving, star-forming universe, we could 
produce precise merger rate predictions without ambi- 
guity. Unfortunately, the dependence dlog Rd / dlog Z 
of rate with metallicity changes dramatically between 
equally plausible models. Extremely strong observa- 
tional constraints are required to limit attention to a 
small region in each parameter and therefore isolate a 
unique Z dependence; e.g., in PSC2 a factor x reduction 
in the parameter volume reduces uncertainty in each pa- 
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rameter by ~ x^^'' . Furthermore, because many of the 
parameters fitted through the comparison to Star Track 
very plausibly could depend implicitly on metallicity, 
such as the strength of stellar winds, a set of parame- 
ters that reproduce the Milky Way need not reproduce 
other star-forming conditions. Thus strong Milky Way 
constraints could but need not eliminate ambiguities as- 
sociated with heterogeneity. 

Strong influences at low metallicity: In the above esti- 
mate we conservatively limit attention to existing pop- 
ulations and employ a fairly narrow metallicity distri- 
bution. Even in the local universe, very young star- 
forming regions can have dramatically lower metallici- 
ties and therefore contribute dramatically more merg- 
ers than allowed for above. Despite their rarity, they 
could dominate the merger rate. Observations of the 
high mass black hole in IC X-10 supports the contention 
that low-metallicity environments of the sort rarely con- 
sidered previously could vastly dominate the present-day 
merger rate. 

5. EXAMPLE: MULTICOMPONENT PREDICTION 

In the above we have argued that an ensemble of bi- 
nary evolution simulations may be needed to generate 
predictions for the distribution of star-forming condi- 
tions within the reach of future gravitational-wave de- 
tectors. A forthcoming paper by Belczynski et al will 
attempt to generate this ensemble in more detail, ex- 
ploring the implications of many different metallicities, 
initial mass functions, and assumptions for binary evolu- 
tion. However, to provide a concrete example that illus- 
trates the challenges associated with heterogeneity, we 
construct merger rate and light predictions for a sim- 
ple two-component universe following the constructive 
procedure in PS-GRB and PS-E. As a sufficiently re- 
alistic example involving an ensemble of metallicities is 
beyond the scope of this paper, we simply adopt choices 
for the metallicity, IMF, and binary evolution model that 
permit us to assemble our illustration from archival cal- 
culations of single-star spectral synthesis and massive 
binary evolution. Specifically, we assume our "ellipti- 
cal" component has low metallicity Z = 0.008 and an 
IMF that at high masses has the fairly flat power law 
din N/ din M = p = —2.125; our "spiral" component 
will have solar metallicity Z = 0.02 and a much steeper 
high- mass power law p = —2.7 (sec iKroupa fc Weidiie^ 
l2003l for an explanation of this choice). The luminos- 
ity density with time is calculated according to [J3] using 
the archived "simple stell ar population" (SSP) models of 
iBruzual fc Charlotl (|2003f ). The merger rate density with 
time is calculated following PS-GRB, adopting random 
but identical assumptions about binary evolution param- 
eters to adopt in the StarTrack model. [Though these 
assumptions are implausible ~ this model assumes much 
higher supernova kicks (cr ^ 950 — 1000 km/s) than are 
currently considered plausible - these models not only 
conveniently involve a m etallicity that appears in the 
IBruzual fc Charlotl (|2003f l archives but also possess ped- 
agogically helpful merger rate histories, as seen below.] 
Finally, following PS-GR B we adopt the two-com ponent 
star formation history of iNagamine "eTaLl (|2f)0l . 

Figure |4] summarizes the results of this concrete ex- 
ample. First, as emphasized in ^ the luminosity and 
merger rate versus time are not simply proportional over- 



all, both because two distinct components (ellipticals and 
spirals) form stars and because light and mergers each 
lag star formation uniquely. Second, depending on the 
type of merging binary of interest, different star form- 
ing conditions can dominate the merger rate. In the 
figure shown, spiral galaxies always dominate the BH- 
NS merger rate; elliptical galaxies dominate the NS-NS 
merger rate; and merging BH-BH binaries are produced 
predominantly in ellipticals early and spirals late. The 
unique response of stars formed in each of the two en- 
vironments, combined the different time-dependent star 
formation histories in each environment, can produce 
many outcomes. Third and not indicated on the figure, 
the characteristic masses of merging binaries generally 
differs in the two components. In the case shown, the 

average detection- weighted chirp mass (M-c / of 

merging BH-BH binaries in ellipticals is similar to that 
of spirals {A.2Mq in spirals, versus 5Mq in ellipticals). 
On the other hand, elliptical galaxies contain noticeably 
less massive merging BH-NS binaries than their spiral 
counterparts (1.9Mq versus 2.6Mq). 

6. CONCLUSIONS 

In anticipation of an era of frequent binary coalescence 
detection and with the goal of divining the limiting as- 
trophysical measurement uncertainties for future obser- 
vations, in this paper we have examined the relevant 
systematic errors intrinsic to proposed absolute normal- 
izations against which gravitational wave detections and 
upper limits can be compared. In other words, we have 
examined the challenges associated with comparing just 
the number of binary merger detections with predictions. 
We find that after a surprisingly small number of detec- 
tions, either much more sophisticated models or richer 
data products (e.g., the observed mass distribution) will 
be needed to further constrain binary evolution. 

For the nearby universe, relevant to initial and en- 
hanced LIGO, we argue that the systematic error asso- 
ciated with using catalog-based normalizations has been 
understated. Though the nominal accuracy of tracers 
of star formation inside a volume, such as blue light as 
adopted in LGC, can be comparable to the systematic er- 
ror target in LIGO (15%), the relevant systematic error 
by adopting a normalization that does not trace old mass 
and remains the same for all binary types and evolution 
models - the error introduced into any comparison be- 
tween the number of detections and predictions - will be 
considerably larger (~ 40%). This systematic error can 
be ameliorated but not eliminated by employing model- 
dependent normalizations. To provide a framework with 
which to calculate this two-band normalization, we in- 
troduced a multi-band galaxy catalog that extends the 
blu e- light catalog presented in LGC; see Figures [1] and 
IA6I We recommend that this catalog and approach be 
applied to re-evaluate the astrophysical systematic errors 
relevant to initial and enhanced LIGO upper limits. 

Advanced detectors will probe the distant universe, for 
which a catalog is impractical. Though merger rates can 
be compared against the average properties of the uni- 
verse, we have demonstrated that treating the universe 
as homogeneous will introduce at least a 40% system- 
atic error, because regions of different metallicity will 
have different relative probabilities of producing massive 
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Fig. 4. — For a two-component universe with the spiral and 
elliptical star formation histories shown in the top panel, predic- 
tions for the time-dependent luminosity density (center panel: U 
(red), B (blue), V (black)) and merg er rate density (bottom panel) 
based on lBruzual fc CharloU I I2003I) and a pair of Star Track pop- 
ulation synthesis models for binary evolution in spiral and ellip- 
tical galaxies that adopt different IMFs and metallicity but oth- 
erwise identical parameters; see the text for details. As in Fig- 
ure |2] merger rates are plotted versus time (t = at present), 
where the peak at t = — 12 Gyr corresponds to a redshift z ~ 2. 
Note the parameters adopted were chosen for convenience in illus- 
tration, not verisimilitude; for example, the very large supernova 
kicks assumed in this mod el are not consisten t with observations 
of isolated pulsars (see.e.e. IHobbs et al.|[2005l) . In each panel the 
contribution overall (solid), from spirals alone (dotted), and from 
ellipticals (dashed) is shown. In the bottom panel, merger rates 
of double neutron star (blue), double black hole (red), and black 
hole neutron star binaries (green) are shown. Note that merger 
rate densities versus time can but need not resemble light versus 
time and that both elliptical and spiral populations can dominate 
a merger rate. 



merging binaries. Wc emphasize our estimate is conser- 
vative, assuming that the only variable in star forma- 
tion is metallicity (e.g., no top-heavy IMFs or alternate 
modes of star formation) and that the universe was al- 
ways homogeneous with a similar metallicity distribution 
to that observed at present. Because binary black hole 
detection rates in particular can be strongly influenced by 
metallicity variations (e.g., due to changes in the initial 
star-final black hole mass relation with metallicity) and 
because black holes arc far more likely to be produced in 
the early universe in the epoch of peak star formation in 
massive galaxies undergoing rapid metallicity evolution 
(binary merger delays for black holes are almost always 
long; based on results in lO'Shaughnessv et al.ll2008al . 
the median merger delay for merging BH-BH binaries 
given steady-state star formation is tbbh — 1 — 3 Gyr. 
depending on assumptions, while for NS-NS binaries it 
is almost always much smaller, tbns ^ 0.3 Gyr), our 
estimate could significantly understate the relevant sys- 
tematic uncertainty. 

To summarize, we recommend the following: (I) When 
interpreting LIGO data as constraints on merger rates, 
unless composition distributions are explicitly incorpo- 
rated into the predictive models, an additional system- 
atic error of order 40% should be included to allow for 
fluctuations in composition and age between galaxies; for 
example, this revised uncertainty will be used in PS-E2 
to explore how advanced LIGO detections might con- 
strain binary merger models. (II) Future merger rate 
predictions should include metallicity evolution and dis- 
tributions, to determine the most likely LIGO detection 
rates when low-metallicity environments are included. 
For example, the forthcoming paper by Belczynski et 
al. will explore evolutionary scenarios over ensemble of 
metallicities in more detail. (Ill) To better assess all 
relevant systematic errors limiting comparisons between 
models and theory, more observational and theoretical 
work is needed to constrain the distribution of fluctua- 
tions, particularly IMF fluctuations early in the universe 
or in clustered star formation. (IV) Finally, to provide 
another handle with which to constrain binary evolu- 
tion, future model constraint papers should describe how 
to compare the detected mass distribution with highly 
model-dependent predictions. Given the immense com- 
putational requirements needed to both thoroughly and 
accurately explore the space of binary evolution mod- 
els, let alone globular clusters, and the relatively mod- 
est benefits that Moore's Law provides to a monte carlo 
simulation sampling a high dimensional space, a care- 
ful balance must be struck between accurately pinning 
down predictions for each model and thoroughly explor- 
ing the model space. The parameter-dependent detection 
efficiency e(D,mi,77i2) and parameter-measurement- 
ambiguity functioi is provided in the gravit ational-wave 
literature (see, e.g. ICutler fc Flanaganlll994[ ) will tell us 
how much we can learn about parameter distributions 
from LIGO and therefore determine where that balance 
will be struck. 

When estimating systematic errors introduced by 
treating the universe as homogeneous, we have for sim- 
plicity assumed all mergers are produced only through 
binary evolution. Interactions in globular clusters are 
expected to be an equally critical channel for form- 
ing merging double black hole binaries; see for example 
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ISadowski et al.l l)2008f ) and references therein. Though 
we have not performed a thorough exploration of pa- 
rameter space, as we were able to do for binary evo- 
lution with StarTrack, wc expect this channel will be 
at least as sensitive to inhomogeneities as isolated bi- 
nary evolution. More critically, this competing channel 
may produce mergers that are indistinguishable from bi- 
nary evolution. The existence of such an unconstrained 
and often indistinguishable channel introduces yet an- 
other large systematic error into interpretation of binary 
compact object detection rates. Further study is critical, 
to determine not only the range of rates these models 
produce in a realistically heterogeneous universe but also 
methods with which to distinguish the randomly-oriented 
and cqual-mass-biascd mergers expected from this chan- 



nel from the more aligned mergers expected from binary 
evolution. 
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Fig. A5. — Lag between light and SFR: Total U (dotted), B ( blue), and V (solid) band l uminosity per Mq of initial star forming mass for 
a starburst at time t = 0, drawn from the Bruzual and Chariot IBruzual fc CharlotI II2003I) spectral synthesis libraries for solar metallicity. 
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APPENDIX 

PHOTOMETRY OF GALAXIES AND BACKGROUNDS 

Photometric conventions: The cumulative luminosities provided in the paper are calculated from the apparent mag- 
nitudes TO and distances d using the solar zero point: 

ix=i0,xlO"°-''(™-*^0'^-^'°8^^''/P")) (Al) 

Specifically, in this paper we adopt for blue light Lq^b ~ 4.7 x 10"^^ erg/s and Mb,q = 5.48; for V-band light 

Q = 4.4 X 10^2 and Mv,q = 4.82; and Lu,q = 1.7 x lO^^ erg/s and Mu,q = 5.66. 
Data sources for catalog: As in LGC, we use a combination of the LED A and TuUy galaxy catalogs to provide corrected 
distances and apparent magnitudes. 

Photometric predictions from SSPs: Rather th an use proportion ality constants that relate the mean star formation 
rate to the present-day light distributions as in iKennicuttI (Il998l). to allow for a more generic comparison we use the 
raw simple stellar population results provided in lBruzual fc CharlotI (|2003[ ) for the kernels Ku{t, Z), Ks ji, Z),Kv{t, Z) 
that relate the star formation rate to the present-day U,B, and V luminosity densities; see Figure [A5l 
Mass-to-light and cumulative mass: In i2]we use the band-to-band differences in (i) the cumulative luminosity distri- 
bution and (ii) fraction from ellipticals to argue that weighting galaxies by light doubtless biases us by of order 20%. 
In the local universe, however, galaxies can be classified by morphology, color, or spectral information into groups with 
roughly similar histories and metallicities. Despite the potential advantage obtained by grouping galaxies with similar 
properties, unless that classification groups galaxies into sufficiently fine groups that all galaxies in that group have a 
similar number of present-day mergers per unit mass, the same limitations often apply; see the discussion in |}21and[3] 
for detailed examples. 

Still, in the spirit of we can introduce another observable that is linearly related to the input star formation rate: 
the mass (with kernel Knjit) = !)■ And as with different bands of light, the merger rate will be well-traced by the mass 
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when the merger rate kernel resembles Km- that is, when it decays very slowly with time. As discussed in fj3l fairly 
few binary evolution models will decay that slowly. Nonetheless, "mass normalization" (treating all star formation 
equally) is a meaningful and extremely complementary normalization to "blue light normalization" (emphasizing only 
the most recent SFR). 

Depending on their stellar content, galaxies can have dramatically different stellar mass to light ratios. The literature 
contains several methods to estimate the relative mass content; for comparison, we adopt two methods, based on 
morphology and color: 

• Morphological classification: iLipunov et al.l ()1995l ) previously used the Tully catalog and a three-component 
morphological classification (into elliptical, spiral, and irregular galaxies) to determine the amount of mass inside 
a sphere at a given radius, using the mass-to-light ratios 

M,/Lb(S) = 1OM0/Lq,b (A2) 
A4/LB(5) = 4.5Af0/L0,B (A3) 
MjLBilrr) = 2Mq/Lq^b (A4) 

where any Sc or Irr galaxy is classified as irregular and young. The fraction of the cumulative "mass" distribution 
obtained with this estimate (Figure EH differs to at least 0(10%) from the cumulative blue luminosity. 

• Color-based M / L estimate: Even galaxies of sim ilar morphological type can differ substantially in their mass to 
light ratios ()Marastonlll 998HBen fc de Jong|l2001h. To esti mate the error in the morphologically-based cumulative 



mass estimate described above, we use iBlanton fc Roweij (|2007t) 's Fi gure 18, which shows a relationship between 
B-V color and il/*/Ly in solar units: 

logio M/Lv « 1.44(5 - V)ab ' 0.76 

-0.3(5 - V)\b (A5) 

{B-V)AB^[B~V)vega-(^-ll (A6) 

where in the first line (following their figure) all magnitudes are referred to an AB magnitude system and in the 
second line an explicit conversion between the two magnitude systems is provided, based on their Table I. 

Based on the differences seen between the cumulative luminosity generated with thi s appr oximation and a 
simple morphological classification (Figure [X6|) or on the spread in IBlanton fc RoweisI (|2007f )'s Figure 18, we 
expect 0(10%) model-dependent uncertainty in the cumulative M and in the fraction of mass contributed from 
ellipticals. 

Asymptotic corrected luminosity per volume: Our catalog consists of extinction-corrected [X =)U,B, V, and far infrared 
(FIR) luminosities. Past the Virgo cluster, the cumulative luminosity Lx{< d) inside a sphere of radius d should revert 
to a mean value 

Air 

where is the extinction-corrected mean galactic emission per unit volume. LGC estimated the mean value jg ~ 
1.98 X 1O~^(1O^°L0^b) Mpc~'^ by correcting the expression in IBlanton et al.l l|2003f l for the luminosity density at 
distances to which advanced detectors will be sensitive (z ~ 0.1) by the expected amount of B-band light that should 
be reprocessed to FIR. 
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Fig. A6. — Mass fraction in ell ipticals (blue) versus distance in kpc. Dotted line: Masses are estimated from B-band luminosity and 
the B-band mass to light ratios of lLipunov et al"] I I1995I) . Solid line: Masses are estimated using each gala xy's visible luminosity (Lv')i 
corrected B-V color, and and Eq. IIASI l. an empirical fit to the data presented in lBlanton fc RoweisI II2007I ) . 



